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Abstract 

We study the phenomenology of a supersymmetric left-right model, assuming minimal supergrav- 
ity boundary conditions. Both left-right and (B-L) symmetries are broken at an energy scale close 
to, but significantly below the GUT scale. Neutrino data is explained via a seesaw mechanism. We 
calculate the RGEs for superpotential and soft parameters complete at 2-loop order. At low energies 
lepton flavour violation (LFV) and small, but potentially measurable mass splittings in the charged 
scalar lepton sector appear, due to the RGE running. Different from the supersymmetric "pure 
seesaw" models, both, LFV and slepton mass splittings, occur not only in the left- but also in the 
right slepton sector. Especially, ratios of LFV slepton decays, such as Br(rR — > fJ>Xi) /Br(r£ — > HXi) 
are sensitive to the ratio of (B-L) and left-right symmetry breaking scales. Also the model predicts 
a polarization asymmetry of the outgoing positrons in the decay —> e + 7, A ~ [0, 1], which differs 
from the pure seesaw "prediction" A = 1. Observation of any of these signals allows to distinguish 
this model from any of the three standard, pure (mSugra) seesaw setups. 
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I. INTRODUCTION 



The most popular explanation for the observed smallness of neutrino masses is certainly 
the seesaw mechanism [lH4|. Literally hundreds of theoretical papers based on "the seesaw" 
have been published since the discovery of neutrino oscillations The seesaw can be 
implemented at tree- level in exactly three realizations {(J: exchange of a fermionic singlet, 
a.k.a. the right-handed neutrino (type-I) [H, 0]; of a scalar triplet (type-II) j^-Ql; or of a 
fermionic triplet (type-Ill) In any of these "seesaw mechanisms" neutrino masses are 
given by m v ~ v 2 /A, where v is the Higgs vacuum expectation value (vev) and A the scale 
of the seesaw. For coefficients 0(1) and A ~ (10 14 — 10 15 ) GeV one finds neutrinos with 
sub-eV masses, just as experimental data demands. Unfortunately, attractive as this idea 
might appear from the theoretical point of view, this estimate also implies that "the seesaw" 
will never be directly tested. 

This situation might change slightly, if supersymmetry (SUSY) is found at the LHC, 
essentially because scalar leptons provide potentially additional information about seesaw 
parameters. Assuming SUSY gets broken at a high energy scale, the seesaw parameters leave 
their imprint on the soft parameters in the Renormalization Group Equation (RGE) running. 
Then, at least in principle, indirect tests of the seesaw become possible^. Indeed, this has 
been pointed out already in [gj], where it was shown that lepton flavour violating (LFV) 
off-diagonal mass terms for sleptons are automatically generated in seesaw (type-I), even 
if SUSY breaking is completely flavour blind at the GUT scale as in minimal supergravity 
(mSugrajE 

Motivated by the above arguments, many authors have then studied LFV in SUSY mod- 
els. For the seesaw type-I, low energy LFV decays such as k — > ^7 and U — > 31 j have been 
calculated in jx — e conversion in nuclei has been studied in |l9l . |20| . The type-II 

seesaw has received much less attention, although it has actually fewer free parameters than 
type-I. The latter implies that ratios of LFV decays of leptons can actually be predicted 



as a function of neutrino angles in mSugra, as has been shown in [21|, |22[. Finally, for 
completeness we mention that LFV in SUSY seesaw type-Ill has been studied in (23)]. 

Measurements at colliders, once SUSY is discovered, can provide additional information. 
LFV decays of left sleptons within mSugra have been studied for type-I in 24j and for type- 
II i n | 22l . |25|. Precise mass measurements might also show indirect effects of the seesaw 



26H28|. Most prominently, type-II and type-Ill seesaw contain non-singlet superfields, so 
gauge couplings run differently from pure MSSM. One then expects that sparticle spectra 
show a characteristic "deformation" with respect to mSugra predictions. From different 
combinations of masses one can form "invariants", i.e. numbers which to leading order depend 
only on the seesaw scale 29|, although there are important corrections at 2-loop 22|, |23[, 
which have to be included before any quantitative analysis can be done. Experimentally 
interesting is also that at the LHC the mass splitting between selectrons and smuons may 

1 In the general minimal supersymmetric extension of the standard model (MSSM) all soft terms are free 

parameters, fixed at the electroweak scale and nothing can be learned about the high energy world. 

2 It might be technically more correct to call this setup the "constrained MSSM" (CMSSM). We will stick 

to the terminology mSugra. 



be constrained down to (9(10~ 4 ) for 30 fb^ 1 of integrated luminosity 30|. In mSugra, one 
expects this splitting to be unmeasurably tiny, whereas in mSugra plus seesaw significantly 
different masses can be generated, as has been shown for type-I in |31| . 

Interestingly, in pure seesaw models with flavour blind SUSY boundary conditions all of 
the effects discussed above show up only in the left slepton sector. Naturally one expects 
that in a supersymmetric model with an intermediate left-right symmetric stage, also the 
right sleptons should contain some indirect information about the high energy parameters. 
This simple observation forms the main motivation for the current paper. Before entering 
in the details of our calculation, let us first briefly discuss left-right symmetric models. 

Quite a large number of different left-right (LR) symmetric models have been discussed in 
the literature. Originally LR models were introduced to explain the observed left-handedness 



of the weak interaction as a consequence of symmetry breaking 32|-|34j. However, LR models 
offer other advantages as well. First, the particle content of LR models contains automat- 
ically the right-handed neutrino and thus the ingredients for generating a (type-I) seesaw 
mechanism Second, the gauge group SU (3)^ x SU(2)l x SU(2)p x U(1)b-l is one °f the 
possible chains through which 5*0(10) 35|, |36| can be broken to the standard model gauge 
group 0. In addition, it has been shown that they provide technical solutions to the SUSY 
CP and strong CP problems 37j and they give an understanding of the U(l) charges of the 
standard model fermions. Interesting only for the supersymmetric versions of LR models, 
(B-L) is gauged and thus, potentially, the low energy theory conserves R-parity j^sf. 

This last argument requires possibly some elaboration. R-parity, defined as Rp = 
(_l)3(B-L)+2s (where _g anc i £ stand for baryon and lepton numbers and s for the spin of the 
particle), is imposed in the MSSM to avoid dangerous baryon and lepton number violating 
operators. However, the origin of Rp is not explained within the MSSM. In early LR models 
SU(2)p doublets were used to break the gauge symmetry. The non-supersymmetric model 



proposed in references [33 
Xl = Xl(1,2, 1, 1) and xr 



34j introduced two additional scalar doublets \l and Xr-, where 
XR (1, 1,2,-1) under SU(3) C x SU(2) L x SU(2) R x U(1) B - L . 



Parity conservation implies that both, \l and Xr, are needed. When the neutral compo- 
nent of xr g e ts a vev, (x°p) ^ 0, the gauge symmetry is broken down to the SM gauge 
group. However, x R is odd under U(1) B -l an d thus, in the SUSY versions of this setup, 
Rp is broken at the same time§. A possible solution to this problem is to break the gauge 
symmetry by SU{2) R fields with even charge under U(1)b-l, i- e - by triplets. For a SUSY 
LR model, this was in fact proposed in reference 4o| . where four triplets were added to the 
MSSM spectrum: A(l, 3, 1, 2), A c (l, 1, 3, -2), A(l, 3, 1, -2) and A c (l, 1, 3, 2). Breaking the 
symmetry by the vev of A c produces at the same time a right-handed neutrino mass via the 
operator L C A C L C , leading to a type-I seesaw mechanism. Depending on whether or not A 



gets a vev, also a type-II seesaw can be generated 41 



However, whether R-parity is conserved in this setup is not clear. The reason is that 

3 Breaking the LR symmetry with triplets can generate also a type-II 

4 Not all 50(10) breaking chains contain a seesaw. Neither does SU(5). It is, of course, straightforward to 
add a seesaw to SU(5). 

5 This could be solved by imposing additional discrete symmetries on the model that forbid the dangerous 
$p operators [39j], but this cannot be regarded as automatic R-parity conservation. 



the minimum of the potential might prefer a solution in which also the right-handed scalar 



42 



Later 



43 



neutrino gets a vev, thus breaking Rp, as has been claimed to be the case in 
calculated some 1-loop corrections to the scalar potential, concluding that Rp conserving 
minima can be found. However, this contradicts the earlier claim 



42 



that 4-loop corrections 



can not eliminate the dangerous Iji p minima. Aulakh and collaborators 



44 



45 



on the other 



hand, showed that by the addition of two more triplets, f2(4, 3, 4, 0) and fl c (l, 4, 3, 0), with 
zero lepton number one can achieve LR breaking with conserved Rp guaranteed already at 
tree-level. Lacking a general proof that the model 40| conserves Rp we will follow 44j, |45 
as the setup for our numerical calculations. 

Finally, for completeness we mention the existence of left-right models with R-parity 
violation. For example, if the left-right symmetry is broken with the vevs of right-handed 
sneutrinos R-parity gets broken as well and the resulting phenomenology is totally different, 
as shown in 46, 47 . 



Compared to the long list of papers about indirect tests of the seesaw, surprisingly little 
work on the "low-energy" phenomenology of SUSY LR models has been done. One loop 
RGEs for two left-right SUSY models have been calculated in [48[ . These two models are 
(with one additional singlet): (a) breaking LR by doublets a la [33, S| and (b) by triplets 
following 40[, but no numerical work at all was done in this paper. The possibility that 
right sleptons might have flavour violating decays in the left-right symmetric SUSY model 
of 40| was mentioned in 49j. A systematic study of all the possible signals discussed above 
for the seesaw case is lacking and to our knowledge there is no publication of any calculation 
of these signals for the model of 44j, |45|. (For completeness we would like to mention that in 
GUTs based on SU(5) one can have the situation the LFV occurs only in the right slepton 

is in a different class from all the 



50 



sector, as pointed out in [501 ] . However, this model 
models discussed above, since it does not contain non-zero neutrino masses.) 

The rest of this paper is organized as follows. In the next section we define the model 



44j . |45| and discuss its particle content and main features at each symmetry breaking scale. 



We have calculated the RGEs for each step complete at the 2-loop level following the general 



description by 5jj using the Mathematica package SARAH 52l454|. A summary is given 



in the appendix, the complete set of equations and the SARAH model files can be found 
Neutrino masses can be fitted to experimental data via a type-I seesaw mechanism 



at 



55 



and we discuss different ways to implement the fit. We then turn to the numerical results. 

for numerical 



56 



The output of SARAH has been passed to the program package SPheno 
evaluation. We calculate the SUSY spectra and LFV slepton decays, such as f^m 
and Tl/r —> ex? and X2 ~^ e A*Xi\ as well as low-energy decays U — > Ijj for some sample 
points as a function of the LR and (B-L) scales. Potentially measurable signals are found 
in both, left and right slepton sectors, if (a) the seesaw scale is above (very roughly) 40 13 
GeV and (b) if the scale of LR breaking is significantly below the GUT scale. Since we find 
sizable LFV soft masses in both slepton sectors, also the polarization in fi — > ey is different 
from the pure seesaw expectation. We then close with a short summary and outlook. 
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TABLE I. Matter content between the GUT scale and the SU(2)r breaking scale. 



II. LEFT-RIGHT SUPERSYMMETRIC MODEL 

In this section we define the model, its particle content and give a description of the 
different symmetry breaking steps. The fit to neutrino masses and its connection to LFV 
violation in the slepton sector is discussed in some detail, to prepare for the numerical results 
given in the next section. We summarize briefly the free parameters of the theory. 



The model essentially follows [44L 145] . We have not attempted to find a GUT completion. 
We will, however, assume that gauge couplings and soft SUSY parameters can be unified, 
i.e. implicitly assume that such a GUT model can indeed be constructed. 



A. Step 1: From GUT scale to SU{2)r breaking scale 

Just below the GUT scale the gauge group of the model is SU(3) C x SU(2)l x SU(2)r x 
U(1)b-l- In addition it is assumed that parity is conserved, see below. The matter content 
of the model is given in table HI Here Q, Q c , L and L c are the quark and lepton superfields 
of the MSSM with the addition of (three) right-handed neutrino(s) v c . 

Two $ superfields, bidoublets under SU(2)l x SU(2)r, are introduced. They contain 
the standard Hd and H u MSSM Higgs doublets. In this model, two copies are needed for a 
non-trivial CKM matrix. Although there are known attempts to build a realistic LR model 



with only one bidoublet generating the quark mixing angles at the loop level 57|, we will 
not rely on such a mechanism. Finally, the rest of the superfields in table [J are introduced 
to break the LR symmetry, as explained above. 

Table [J shows also the gauge charges for the matter content in the model. In particular, 
the last column shows the B — L value for the different superfields. However, the following 
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definition for the electric charge operator will be used throughout this paper 

Q = ISL + hR + (1) 



and thus the U{1)b-l charge is actually ^y ± - 

With the representations in tabled the most general superpotential compatible with the 
gauge symmetry and parity is 

W = Y Q Q$Q C + Y L LQL C - + fLAL + f*L c A c L c 
+ aAQA + a*A c Q c A c + + a*fi c $$ 

+ Ma A A + M A A C A C + M n ttn + M^fL c VL c . (2) 

Note that this superpotential is invariant under the parity transformations Q -H- (Q c )*, 
L <-> (L c )*, $ <-> A ^ (A c )*, A ^ (A c )*, ft (fi c )*. This discrete symmetry fixes, for 
example, the L C A C L C coupling to be /*, the complex conjugate of the LAL coupling, thus 
reducing the number of free parameters of the model. 

Family a nd g auge indices have been omitted in eq. (121), more detailed expressions can 
be found in 44j. Yq and Yl are quark and lepton Yukawa couplings. However, with two 
bidoublets there are two copies of them, and thus there are four 3x3 Yukawa matrices. 
Conservation of parity implies that they must be hermitian. \i is a 2 x 2 symmetric matrix, 
whose entries have dimensions of mass, / is a 3 x 3 (dimensionless) complex symmetric 
matrix, and a is a 2 x 2 antisymmetric matrix, and thus it only contains one (dimensionless) 
complex parameter, ol\i. The mass parameters Mq and Ma can be exchanged for vr and 
Vbli the vacuum expectation values of the scalar fields that break the LR symmetry, see 
below. 

The soft terms of the model are 

- C soft = m 2 Q Q j Q + m 2 Qc Q cj Q c + m\l)L + m\ c L c] L c 

+ m|$ f $ + m^A f A + m\A ] A + m\ c A c] A c + m| c A ct A c 

+ ml tfn + m 2 n c tt c1 tt c + 1 [M^B + M 2 {W L W L + W R W R ) + M 3 gg + h.c] 

+ [T Q Q$Q C + T L L$L C + TfLAL + T* f L c A c L c 

+ T a AttA + T*A c Vl c A c + T Q 1M>$ + T*fi c $$ + h.c] 

+ [B^ + B Ma AA + B Ma *A c A c + B Ma Sm + 5 Mn *fi c fi c + h.c] . (3) 

Again, family and gauge indices have been omitted for the sake of simplicity. The LR 
model itself does not, of course, fix the values of the soft SUSY breaking terms. In the 
numerical evaluation of the RGEs we will resort to mSugra-like boundary conditions, i.e. 
^0^3x3 = m Q = m Qc = m l = m L^ m o I 2x2 = tti%, m\ = m\ = m\ = m\ c = m\ c = ml = 
m 2 nc , Mi /2 = Mi = M 2 = M 3 , T Q = A Y Q ,T L = A Y L ,T f = A f,T a = A a,T a = A a, 
Bp = B ,B Ma = BqM^, B Mn = B Mq. The superpotential couplings /, Yq and Y L are 
fixed by the low-scale fermion masses and mixing angles. Their values at the GUT scale are 
obtained by RGE running. This will be discussed in more detail in section III Dl 



(3 



The breaking of the LR gauge group to the MSSM gauge group takes place in two steps: 
SU(2) R x U(1)b-l — > U(1)r x U(V)b-l — > U(l)y- In the first step the neutral component 
of the triplet Q takes a vev: 

(fi c °) = ^ (4) 

which breaks SU(2) R . However, since I3r(Q c0 ) = there is a U(1) R symmetry left over. 
Next, the group U(1) R x U(1)b-l is broken by 

(A c0 ) = ^ , (A c0 ) = ^ . (5) 



The remaining symmetry is now U(l)y with hypercharge defined as Y = I 3R + 



2 - 



The tadpole equations do not link f2 c , A c and A c with their left-handed counterparts 



441 and the 



due to supersymmetry. Thus, the left-handed triplets can have vanishing vevs 
model produces only a type- 1 seesaw. 

Although a "hierarchy" between the two breaking scales may exist, Vbl *C v r , one cannot 
neglect the effects of the second breaking stage on the first one, since mass terms of Q and 
A enter in both tadpole equations. If we assume Vbl — v bl the tadpole equations of the 
model can be written 

^- = A\M n \ 2 v R + ^\a\ 2 v 2 BL v R - l -v\ L [a* (Ma + M n ) + c.c] = , (6) 

= \M a \ 2 v B l + \\a\ 2 {v 2 BL + v R )v BL - \v BL v R [a* (Ma + M n ) + c.c] = . (7) 

In these equations (small) soft SUSY breaking terms have been neglected. Similarly, at 
this stage there are no electroweak symmetry breaking vevs and v u . From equations 
(jSJ) and one sees that, in fact, there is an inverse hierarchy between the vevs and the 
superpotential masses Ma, Mq, given by 

v R = ^ , ^ = -(2M A M n )^ . (8) 



And so, Vbl ^ v r requires Ma 3> Mq, as has already been discussed in 44j . 



B. Step 2: From SU(2) R breaking scale to U(1)b-l breaking scale 

At this step the gauge group is SU (3) c x SU (2) L xU(l) R xU (1) b-l- The particle content 
of the model from the SU(2) R breaking scale to the U(1)b-l breaking scale is given in table 

El 

Some comments might be in order. Despite Ma being of the order of v R (or larger), see 
eq.flB]), not all components of the A superfields receive large masses. The neutral components 
of A c and A c lie at the v R l scale. One can easily check that the F-term contributions to 
their masses vanish in the minimum of the scalar potential eq. ((§]). Moreover, Q c does not 
generate D-terms contributions to their masses. Therefore, contrary to the other components 
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TABLE II. Matter content from the SU(2)r breaking scale to the U{1)b-l breaking scale. 



of the A triplets, they only get masses at the Vbl scale. On the other hand, one might guess 
that all components in the Q,Q C superfields should be retained at this stage, since their 
superpotential mass is required to be below Vbl- However, some of their components get 
contributions from SU (2)# breaking, and thus they become heavy. The charged components 
of Q c do develop large masses, in the case of the scalars through D-terms, while in the case 
of the fermions due to their mixing with the charged gauginos W R , which have masses 
proportional to vr. However, the neutral components of Q c do not get SU(2)r breaking 
contributions, since they have I^r{VL cQ ) = 0, and then they must be included in this energy 



regime. See reference |45| for a more quantitative discussion. 

After SU(2)r breaking the two bidoublets $i and $2 get split into four SU(2)l dou- 
blets. Two of them must remain light, identified with the two Higgs doublets of the MSSM, 
responsible for EW symmetry breaking, while, at the same time, the other two get masses 
of the order of vr. This strong hierarchy can be only obtained by imposing a fine-tuning 
condition on the parameters involved in the bidoublet sector. 

The superpotential terms mixing the four SU(2)l doublets can be rewritten as 

W M = (H f d ) T M H Hl (9) 

where H d = (H},H%) and H[ = (H^H^) are the flavour eigenstates. In this basis reads 
the matrix 

\ H12 ~ a 12 M R fi 22 J 
where the relations = fiji and = — <x,-j have been used and Mr = ^ has been defined. 
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In order to get two light doublets we impose the fine-tuning condition 



45 



Det(Mff) = nuim - (pis - ai 2 M 2 R ) = . (11) 

The result of eq. (ITTj) is to split the two Higgs bidoublets into two pairs of doublets (Hd, H u )l 
and (Hd, H U ) R , where (Hd, H u )l is the light pair that appears in table [TTJ and (Hd, H U ) R 
a heavy pair with mass of order of v R . In practice, equation (fTTj) implies that one of the 
superpotential parameters must be chosen in terms of the others. Since this fine-tuning 
condition is not protected by any symmetry, the RGEs do not preserve it, and one must 
impose it at the SU(2) R breaking scale. In our computation we chose to compute in 
terms of the free parameters \i\ 2 , &12 and vr. 

In order to compute the resulting couplings for the light Higgs doublets one must rotate 
the original fields into their mass basis. Since Mh is not a symmetric matrix (unless a.12 = 0) 
one has to rotate independently H d and H£, i.e. H d = DHJ 1 , H[ = UH™, where D and U 
are orthogonal matrices and H™ = (H d , Hf) and H™ = (H^, H^) are the mass eigenstates. 
This way one finds 

W M = {H f d ) T M H Hl = (H?) T D T M H UH™ = (H™) T M H H™ (12) 
where Mjj is a diagonal matrix, with eigenvalues 
M 2 H>1 = , 

M 2 H ,2 = 4" + 2a 2 12 M 2 R (A ~ /4) + (A + £2?) ■ (13) 

P22 

The D and U rotations are, in general, different and we parametrize them as 

D _( cos6>! sin^i \ U—( cos ^ 2 sin ^ 2 | ^4) 

I — sin 9\ cos 0\ J ' l— sin 9 2 cos 62 J 



and get 



H] = cos 0iH% + sin 6 1 H d 

r2 _ , ■ , n ttL 1 n tjR 



R 

d 

H'i = - sin d x H% + cos Q X H% , (15) 



and similar for H u . In general the angles Q\ and 9 2 are different. However, they are connected 
to the same matrix Mh and can be calculated by diagonalizing M#(M#) T or (M#) T M# 
and one finds 

tan^ = ^*. (16) 

A*22 

In these expressions Det(M#) = has been used to simplify the result. Exact Det(M#) = 
implies that the /i-term of the MSSM is zero, so this condition can only be true up to small 
corrections, see the discussion below. Note that there are two interesting limits. First, 
fii 2 ^> ai 2 M R : this implies tan#i = tan 9 2 and therefore D = U. This is as expected, since 
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that limit makes Mh symmetric. And, second, [i\ 2 <C ol\ 2 Mr : this implies tan#i = — tan 9 2 
and therefore D = U T . 

The superpotential at this stage is 

W = Y u QH u u c + Y d QH d d c + Y e LH d e c + Y u LH u u c + fiH u H d 
+ / c Vz/A c0 + M^ c A c0 A c0 + aA^A + a^A c0 A c0 fi c0 
+ bttH d H u + b c tt c0 H d H u + M n nn + M n cfi c0 fi c0 . (17) 

Particles belonging to the same SU(2)r gauge multiplets split due to their different U{1)r 
charges. At this stage both the LR group, that symmetrizes the SU{2)l and SU(2)r gauge 
interactions, and the discrete parity symmetry that we imposed on the couplings are broken. 
The soft terms are 

2 i7cK7c _i_ tti 2 Ac* Ac 1 m 2 f t f — 1 — ™2 ~ c t ~ c — 1 — m 2 ~ c t ~, 



>C so /t = m 2 nQ^Q + m 2 uC u c u c + m^ c <i c <i c + m 2 L lJ L + m^e c e c + m^. 



+ m\E\R u + m 2 Hd H d H d + m^ A c0T A c0 + m| c0 A c °tA c0 

+ mffin + m 2 nc0 n c01 n c0 + - [M^B + M L W L W L + M R W° R W° R + M 3 ^ + /i.e.] 

+ [T u QH u u c + T d QH d d c + T e LH d e- + T u LH u v° (18) 
+ T} c u^A c0 + T* c A c0 fi c0 A c0 + T b QH d H u + T b c{l c0 H d H u + /i.e.] 
+ [£^if u # d + 5 M i c A c0 A c0 + 5 A/f2 fifi + B M oQ c0 n c0 + h.c] . 

Again we suppress gauge and family indices. 

We must impose matching conditions at the SU{2)r breaking scale. These are for super- 
potential parameters given by 

Y d = Y^ cos 9 1 - Yq sin 9 1 , Y u = —Yq cos 9 2 + Y 2 sin 9 2 , 
Y e = Y[ cos 9 l - Yl sin 9 1 , Y v = -Y[ cos 9 2 + Y 2 sin 9 2 , 



at 



Jc - J , -c ^ 

Mi c = Ml , M^ c = M* , 

b = 2aR, b c = V2a*R, (19) 
where i? = sin(6' 1 — 9 2 ). For the soft masses we have 

m lc = ™% = m Q- , (20) 
m 2 eC = m 2 vC = m 2 Lc , 
2 2 

TYlj^cO = 7Tlj\c , 

2 2 

= ITT'^c ; 

2 _ 2 

m Qc0 — m Qc , 
M L = M R = M 2 . 

Soft trilinears matching follow corresponding conditions. In addition, one has 

m 2 „ = cos^K),, + sin^.K)^ - sinftcosO, [(m|) 12 + (m|) 21 ] , 
< = cos 2 * 2 «)n + fcW)- - si. « 2 c OS 2 [«) 12 + «) 21 ] , 
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as obtained when the operator m|$t$ is projected into the light Higgs doublets operators 
(HfflH% and {H^yH^. Gauge couplings are matched as <?l = 9r = 92- 



C. Step 3: Prom U(l) B -L breaking scale to EW/SUSY scale 

We mention this stage only for completeness, since the last regime is just the usual MSSM. 
We need matching conditions in the gauge sector. Since U(1)r x U(1)b-l breaks to U(1)y, 
the MSSM gauge coupling g\ will be a combination of g^ and gsL- The resulting relationship 
is 

V59R9BL / 01 v 

9i = 7= 2 = 2 : • 21 
V 2 9r + 39bl 

Analogously, the following condition holds for gaugino masses 



M^MSSM) = 2q2rM \ + 3g f Mr . (22) 
2 9r + 3g BL 



Note that in the last two equations the gauge couplings are GUT-normalized. Electroweak 
symmetry breaking occurs as in the MSSM. We take the Higgs doublet vevs 

(H° d ) = ^ , (#°> = ^ , (23) 

as free parameters and then solve the tadpole equations to find /xmssm and B* 1 . /xmssm must 
be different from zero, that is Det(M#) can not be exactly zero. Instead the tuning must be 
exact up to Det(M#) = CK/Wmssm)- usu& 1 tan/3 = ^ is used as a free parameter. Also 
the sign of /xmssm is not constrained as usual. 



D. Neutrino masses, LFV and Yukawa couplings 

Neutrino masses are generated after U(1)b-l breaking through a type-I seesaw mecha- 
nism. The matrix /* leads to Majorana masses for the right-handed neutrinos once A c0 
gets a vev. We define the seesaw scale as the lightest eigenvalue of 

M s = f l c VBL ■ (24) 

As usual, we can always rotate the fields to a basis where Ms is diagonal. However, 
this will introduce lepton flavour violating entries in the Y^. Yukawas, see discussion below. 
As mentioned above, contrary to non-supersymmetric LR models |2[, there is no type-II 
contribution to neutrino masses. 

Global fits to all available experimental data provide values for the parameters involved 



in neutrino oscillations, see table II III for updated results and ref. 59|, |60| for experimental 



results. As first observed in 



61 



these data imply that the neutrino mass matrix can be 
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parameter 


best fit 


2-0- 


A 772,91 flO- 5 eV 2 l 
|Am§ 1 |[l(r 3 eV 2 ] 

sin 2 612 

sin 2 6*23 

sin 2 6»i3 


— U.lo 

2 4n+ - 12 
n Qi 0+0.019 

u - ou -0.06 
U.UlO_ 009 


7.22 - 8.03 
2.18 - 2.64 
0.29 - 0.36 
0.39 - 0.63 
< 0.039 



TABLE III. Best-fit values with l-cr errors and 2-cr intervals (1 d.o.f.) taken from the reference [58| . 
which is updated continuously on the web. 



diagonalized to a good approximation by the so-called tri-bimaximal mixing pattern: 



U' 



TBM 



( 


fi 


fl 





\ 




V 3 


V 3 






1 


1 


1 






V6 




V2 




K 


1 


1 


1 


) 


V6 




V2 



(25) 



The matrix product Y v ■ {fl)' 1 ■ Yj is constrained by this particular structure. LFV entries 
can be present in both Y v and fl, see also the discussion about parameter counting in the 
next subsection. However, in the numerical section we will consider only two specific kinds 
of fits: 



• Y^-fit: flavour structure in Y v and diagonal j\. 

• /-fit: flavour structure in fl and diagonal Y v . 

While at first it may seem either way of doing the fit is equivalent, fl and Y v in our setup can 
leave different traces in the soft slepton mass parameters if vbl vr. This last condition 
is essential to distinguish between both possibilities, because otherwise one obtains the 
straightforward prediction that LFV entries in left and right slepton are equal, due to the 
assumed LR symmetry above Vr. 



These two types of fit were already discussed in reference 62j, which investigates low 
energy LFV signatures in a supersymmetric seesaw model where the right-handed neutrino 
mass is generated from a term of the form fA c u c u c . When the scalar component of A c 
acquires a vev a type-I seesaw is obtained, generating masses for the light neutrinos. There- 
fore, this model has the ingredients to accommodate a Y^-fit, named as Dirac LFV in j^ . 
or a /-fit, named as Majorana LFV. Note, however, that the left-right symmetry, central in 
our work, is missing in this reference, thus implying different signatures at the electroweak 
scale. 

The difference in phenomenology of the two fits can be easily understood considering 
approximated expressions for the RGEs for m 2 L and m 2 c . In the first step, from the GUT 
scale to the v R scale RGEs at 1-loop order can be written in leading-log approximation as 
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Ami = — (3//* + vfyf (3 ml + K) In (^f ) , 

Am|. = (3/V + V<") (3ml + 4) In (^) . (26) 

Of course, also the A parameters develop LFV off-diagonals in the running. We do not give 
the corresponding approximated equations for brevity. After parity breaking at the Vr scale 
the Yukawa coupling Yi splits into Y e , the charged lepton Yukawa, and Y u , the neutrino 
Yukawa. The later contributes to LFV entries in the running down to the Vbl scale. Thus, 

Ami ~ ~^ Y ^ Kk + Ai\ VR ) In (j^ , 

Am 2 eC ~ , (27) 

where m 2 L \ VR is the matrix m\ at the scale vr and A 2 e \ VR is defined as T e = Y e A e and also has 
to be taken at Vr. In order to understand the main difference between the two fits, let us first 
consider the /-fit. This assumes that Y v is diagonal at the seesaw scale and thus the observed 
low energy mismatch between the neutrino and charged lepton sectors is due to a non-trivial 
flavour structure in j\. Of course, non-diagonal entries in / generate in the running also 
non-diagonal entries in Y v and Y e , but these can be neglected in first approximation. In this 
case, equations f l26j) and (127]) show that the LR symmetry makes m\ and m\ c run with the 
same flavour structure and the magnitudes of their off-diagonal entries at the SUSY scale 
are similar. If, on the other hand, Y v is non-trivial (Y^-fit), while / is diagonal, the running 
from the GUT scale to the vr scale induces again the same off-diagonal entries in m\ and 
m\c. However, from vr to vbl the off-diagonals entries in m\ continue to run, while those 
in m\a do not. This effect, generated by the right-handed neutrinos via the Y v Yukawas, 
induces additional flavour violating effects in the L sector compared to the R sector. Seeing 
LFV in both left and right slepton sectors thus allows us to indirectly learn about the high 
energy theory. We will study this in some detail in the numerical section below. 

E. Parameter counting 

Let us briefly summarize the free parameters of the model. With the assumption of 
mSugra (or better: mSugra-like) boundary conditions, in the SUSY breaking sector we 
only have the standard parameters m , Mu2, A , tan/3, sign^usSM)- Thus, we count 4+1 
parameters in the soft terms. We note in passing that the soft terms of the heavy sector, 
of course, do not have to follow strictly the conditions outlined in equation ([3]), as long as 
these parameters are small compared to vbl there are no changes compared to the above 
discussion. 

In the superpotential we have a, a, /i, and Mj). This leaves, at first sight, 7 parameters 
free. However, we can reduce them to 4+2 parameters as follows. S only 
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contains one free parameter: a\2- The matrix /i has 3 entries, but one of them, /in, is fixed 
by the fine-tuning condition Det(M^) = (^(/^mssm)- This leaves two free parameters, /ii2, 
//22- We have traded Ma and Mc> for the vevs vr, Vbl, since ln(^-) and ln( ^ T ) enter 
into the RGEs and thus can, at least in principle, be determined from low-energy spectra. 
There are then in summary 6 parameters, four independent of low-energy constraints and 
two which could be fixed from LFV data, see below. 

In addition, in the superpotential we have the Yukawa matrices Yq., Y^. and /. Let's 
consider the quark sector first. Since we can always go to a basis in which one of the Yq. is 
diagonal with only real entries, there are 12 parameters. Ten of them are fixed by six quark 
masses, three CKM angles and the CKM phase, leaving two phases undetermined. 

In the lepton sector we have the symmetric matrices, and Yl 2 . As with the quark 
sector, a basis change shows that there are only 12 free parameters. / is symmetric and thus 
counts as another 9 parameters. Going to a basis in which / is diagonal does not reduce the 
number of free parameters, since in this basis we can no longer assume one of the Y^. to be 
diagonal. In summary there are thus free 21 parameters in these three matrices. 

In the simple, pure seesaw type-I with three generations of right-handed neutrinos the 
number of free parameters is 21. Only 12 of them can be fixed from low-energy data: three 
neutrino and three charged lepton masses, three leptonic mixing angles and three phases (two 
Majorana and one Dirac phase). However, as pointed out in ll|, in principle, m\ contains 
9 observable entries and thus, if the normalization (i.e. m , A , tan/3 etc.) is known from 
other sfermion measurements, one could re-construct the type-I seesaw parameters §. 

How does the SUSY LR model compare to this? We have, as discussed above, also 
21 parameters in the three coupling matrices, but neutrino masses depend also on Vbl- 
However, in principle, we have 9 more observables in m^, assuming again that the soft 
SUSY breaking terms can be extracted from other measurements. Since in the RGEs also 
Vr appears we have in total 23 parameters which need to be determined. The number of 
observables, on the other hand is fixed to 30 in total, as we have 12 (low-energy lepton 
sector) plus 9 (left sleptons) plus 9 (right sleptons) possible measurements. 



III. NUMERICAL RESULTS 
A. Procedure for numerics 

All necessary, analytical expressions were calculated with SARAH. For this purpose, two 
different model files for the model above the two threshold scales were created and used 
to calculate the full set of 2-loop RGEs. SARAH calculates the RGEs using the generic 
expressions of in the most general form respecting the complete flavour structure. These 
RGEs were afterwards exported to Fortran code and implemented in SPheno. As starting 
point for the RGE running, the gauge and Yukawa couplings at the electroweak scale are 

6 Of course, this discussion is slightly academic, since at least one of the Majorana phases will never be 
measured in praxis. 
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used. In the calculation of the gauge and Yukawa couplings we follow closely the procedure 
described in ref. 



56 



the values for the Yukawa couplings giving mass to the SM fermions 
and the gauge couplings are determined at the scale Mz based on the measured values for 
the quark, lepton and vector boson masses as well as for the gauge couplings. Here, we 
have included the 1-loop corrections to the mass of W- and Z-boson as well as the SUSY 
contributions to Svb f° r calculating the gauge couplings. Similarly, we have included the 



complete 1-loop corrections to the self-energies of SM fermions 
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Moreover, we have 



resummed the tan /3 enhanced terms for the calculation of the Yukawa couplings of the b- 



quark and the r-lepton as in |56|. The vacuum expectation values v& and v u are calculated 
with respect to the given value of tan/3 at Mz- Since we are working with two distinct 
threshold scales, not all heavy fields are integrated out at their mass and the corresponding 1- 
loop boundary conditions at the threshold scales are needed. It is known that these particles 
cause a finite shift in the gauge couplings and gaugino masses. The general expressions are 
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^(^i^^KS))' (28) 
M ^ M <( i± ^w k (3)) ■ (29) 

I\{r) is the Dynkin index of a field transforming as representation r with respect to the 
gauge group belonging to the gauge coupling g^ M is the mass of this particle and Mt is the 
threshold scale. When evaluating the RGEs from the low to the high scale, the contribution 
is positive, when running down, it is negative. The different masses used for calculating 
the finite shifts are the eigenvalues of the full tree-level mass matrix of the charged, heavy 
particles removed from the spectrum. The correct mass spectrum is calculated in an iterative 
way. The GUT scale is defined as the scale at which gsL = 92 = 9gut holds. Generally, 
there is difference with g^ to gcuT in the percent range, the actual numerical mismatch 
depending on the scales Vbl and Vr and being larger for lower values of Vbl and vr. It 



has been stressed in particular in |65| that within supersymmetric LR models, the LR 
symmetry breaking scale has to be close to the GUT scale, otherwise this mismatch will 
grow too large. However, in j66| it was pointed out that, among other possibilities, GUT 
thresholds - unknown unless the GUT model, including the complete Higgs sector used to 
break the GUT symmetry, is specified - can lead to important corrections, accounting for 
this apparent non-unification. (For a discussion of these effects in the context of SU(5) see 



671].) We simply use g^L — 92 = 9gut and attribute departures from complete unification 
to (unknown) thresholds. After applying the GUT scale boundary conditions, the RGEs are 
evaluated down to the low scale and the mass spectrum of the MSSM is calculated. The 
MSSM masses are, in general, calculated at the 1-loop level in the DR scheme using on-shell 
external momenta. For the Higgs fields also the most important 2-loop contributions are 
taken into account. We note that the corresponding Fortran routines are also written by 
SARAH but they are equivalent to the routines included in the public version of SPheno 
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FIG. 1. Example of spectra at the SUSY scale and its dependence on vrl (left side) and vr (right 
side). The masses of four states are shown: x.i (blue line), X2 (blue dashed line), /Ir (red line) 
and /ij, (red dashed line). In both panels the mSugra parameters have been taken as in the SPS3 
benchmark point. 

based on j63|. The iteration stops when the largest change in the calculation of the SUSY 
and Higgs boson masses at Msusy is below one per-mille between two iterations. 



B. Mass spectrum 



The appearance of charged particles at scales between the electroweak scale and the GUT 



scale leads to changes in the beta functions of the gauge couplings 
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This does not 



only change the evolution of the gauge couplings but also the evolution of the gaugino and 



scalar mass parameters 
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The LR model contains additional triplets, and similar to 



23| the mass spectrum at low energies is shifted with 



what is observed in the seesaw models 
respect to mSugra expectations. Two examples of this behaviour are shown in figure [U In 
this figure we show the two lightest neutralino masses and the masses of the left and right 
smuons versus vrl (left side) and Vr (right side). We note that also all other sfermion and 
gaugino masses show the same dependence and in general smaller values are obtained for 
lower values of vrl and Vr. One finds that gaugino masses depend stronger on Vrl an d Vr 
than sfermion masses and that right sleptons are the sfermions for which the sensitivity to 
these vevs is smallest. 

The change in the low energy spectrum, however, maintains to a good degree the standard 
mSugra expectation for the ratios of gaugino asses, as shown in figures [2] and |3j Here, figure 
|2] shows the ratios M\jM.% and M2/M3 versus vrl, while figure [3] shows the same ratios 
versus vr. Shown are the results for three different SUSY points, which in the limit of 
vr,v bl — > m GUT approach the standard SPS points SPSla' (68|, SPS3 and SPS5 



69 



For 
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FIG. 2. Gaugino mass ratios as a function of vrl for the fixed value vr = 10 GeV. To the left, 
M1/M2, whereas to the right M2/M3. In both figures the three coloured lines correspond to three 
mSugra benchmark points: SPSla' (blue), SPS3 (green) and SPS5 (red). Note the small variation 
in the numbers on the Y axis. 

example, the ratio M\jMi is expected to be (5/3) tan 2 9w — 0.5 at 1-loop order in mSugra. 
The exact ratio, however, depends on higher order corrections, and thus on the SUSY 
spectrum. The LR model will thus appear rather mSugra like, if these ratios are measured. 
Only with very high precision on mass measurements, possible only at a linear collider, can 
one hope to find any (indirect) dependence on vrl and vr. 



C. LFV of leptons 

Lepton flavour violation in charged lepton decays has attracted a lot of attention for 
decades. Processes like \x — > are highly suppressed in the standard model (plus non-zero 



neutrino masses) due to the GIM mechanism |70|, and thus the observation of these rare 
decays would imply new physics. The MEG experiment [ti| is currently the most advanced 
experimental setup in the search for fi + — > e +r y. This rare decay will be observed if its 
branching ratio is above the MEG expected sensitivity, around Br(fi — > e^) ~ 10~ 13 . 

LFV decays like /j — > ^7 are induced by 1-loop diagrams with the exchange of neutralinos 
and sleptons. They can be described by the effective Lagrangian, see for example the review 

777 • — 

£eff = e^ka^{AlP L + AlP R )lj + h.c. . (30) 

Here Pl,r = |(1 =F 7s) are the usual chirality projectors and therefore the couplings Al 
and Ar are generated by loops with left and right sleptons, respectively. In our numerical 
calculation we use exact expressions for Al and Ar. However, for an easier understanding 
of the numerical results, we note that the relation between these couplings and the slepton 
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v R [GeV] 



v R [GeV] 



FIG. 3. Gaugino mass ratios as a function of vr for the fixed value vbl = 10 14 GeV. To the left, 
M1/M2, whereas to the right M2/M3. In both figures the three coloured lines correspond to three 
mSugra benchmark points: SPSla' (blue), SPS3 (green) and SPS5 (red). Note the small variation 
in the numbers on the Y axis. 



soft masses is very approximately given by 



{™L)i 



m 



Imt 



SUSY 



ml 



(31) 

'■SUSY 

where msusY is a typical supersymmetric mass. Here it has been assumed that (a) 
chargino/neutralino masses are similar to slepton masses and (b) A-terms mixing left- 
right transitions are negligible. Therefore, due to the negligible off-diagonal entries in 77ig C , 
a pure seesaw model predicts Ar ~ 0. 

The branching ratio for — > ^7 can be calculated from the previous formulas. The result 

is 

Br(k = -^f- (\A l [\ 2 + \A^\ 2 ) Br(h Ijvft) . (32) 



Figure H] shows two examples for Br(n — > ey) in the tuq, M1/2 plane. Here, we have fixed 
Vbl = 10 14 GeV and v R = 10 15 GeV and show to the left M s = 10 12 GeV, whereas to the 
right Ms = 10 13 GeV. Here we have assumed a degenerate spectrum right-handed neutrinos 
which we denote by Ms = Mm. Once Yukawas are fitted to explain the observed neutrino 
masses, the branching ratio shows an approximately quadratic dependence on the seesaw 
scale, with lower Ms giving smaller Br([i — > &-/). As expected, the branching ratio also 
strongly decreases as mo and/or M1/2 increase. This is because the superparticles in the 
loops leading to fi — > &y become heavier in these directions, suppressing the decay rate. In 
fact, from equations ( IHTj) and ( 132]) one easily finds the dependence 



Br(ji — > ery) 



r<2 



nr 



(33) 



SUSY 
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M Seesaw= 10 ( GeV ) tanp=10, A Q =0 (GeV) 



M Seesaw= ™™ (GeV) tanp=10, A =0 (GeV) 
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FIG. 4. Contours of Br([i -> e-f) in the m ,M 1 / 2 plane for vrl = 10 GeV and vr = 10 15 GeV. 
To the left Ms = 10 12 GeV, whereas to the right M5 = 10 13 GeV. Neutrino oscillation data have 
been fitted with the Y,, fit. 



which shows that Br(fi — > ej) decreases as Tng^ SY . 

It is also remarkable that for a given seesaw scale, Br(fi is sizeably larger in the 

LR model than in a pure seesaw type-I model, see for example |25| . The explanation of this 
is that right sleptons contribute significantly in the LR model to Br(fi — > cy) and these 
contributions are absent in seesaw models. 

As already discussed, a pure seesaw model predicts simply Ar ~ 0. However, in the LR 
model we expect a more complicated picture. Left-right symmetry implies that, above the 
parity breaking scale, non-negligible flavour violating entries are generated in ra^. Therefore, 
Ar 7^ is obtained at low energy. The angular distribution of the outgoing positron at, for 
example, the MEG experiment could be used to discriminate between left- and right-handed 



polarized states 73|, |74J. If MEG is able to measure the positron polarization asymmetry, 
defined as 





2 _ 


Ar 


2 


\A L 


2 + 


Ar 


2 



-+ e+ 7 ) = (34) 



there will be an additional observable to distinguish from minimal seesaw models. In a pure 
seesaw model one expects A — +1 to a very good accuracy. However, the LR model typically 
leads to significant departures from this expectation, giving an interesting signature of the 
high energy restoration of parity. 

Figure [5] shows contours for A(fi + —> e + 7) in the mo, M1/2 plane. For the corresponding 
branching ratios see figure HI Note the rather strong dependence on m . The latter can be 
understood as follows. Since vrl in these examples is one order of magnitude smaller than 
vr, and the Y v fit has been used, the LFV mixing angles in the left slepton sector are larger 
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FIG. 5. Contours of A(fJ, — > e + j) in the mo,Mi/ 2 plane. To the left Ms = 10 GeV, whereas to 
the right Ms = 10 13 GeV. The parameters have been chosen as in figure HI 



than the corresponding LFV entries in the right sleptons. At very large values of tuq, were 
the masses of right and left sleptons are of comparable magnitude, therefore "left" LFV is 
more important and the model approaches the pure seesaw expectation. At smaller values 
of mo, right sleptons are lighter than left sleptons, and due to the strong dependence of 
fi — > e7 on the sfermion masses entering the loop calculation, see eq. ( 13TT) . Ar and Al can 
become comparable, despite the smaller LFV entries in right slepton mass matrices. In the 
limit of very small right slepton masses the model then approaches A ~ 0. We have not 
explicitly searched for regions of parameter space with A < 0, but one expects that negative 
values for A are possible if vbl is not much below Vr and sleptons are light at the same 
time, i.e. small values of mo and M\ji. Note that, again due to the LR symmetry above to 
Vr, the model can never approach the limit A — — 1 exactly. 

The positron polarization asymmetry is very sensitive to the high energy scales. Figure 
[6] shows A as a function of vr for Ms = 10 13 GeV, vbl — 10 14 GeV and the mSugra 
parameters as in the SPS3 benchmark point. The plot has been obtained using the Y v fit. 
This example shows that as Vr approaches rricuT the positron polarization A approaches +1, 
which means Al dominates the calculation. This is because, in the Y v fit, the right-handed 
LFV soft slepton masses, and thus the corresponding Ar coupling, only run from itigut to 
vr. 

A(fi + — > e + 7) also has an important dependence on the seesaw scale. This is shown 
in figure [TJ where A is plotted as a function of the lightest right-handed neutrino mass. 
This dependence can be easily understood from the seesaw formula for neutrino masses. It 
implies that larger Ms requires larger Yukawa parameters in order to fit neutrino masses 
which, in turn, leads to larger flavour violating soft terms due to RGE running. However, 
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FIG. 6. Positron polarization asymmetry A(fi + —¥ e + 7) as a function of vr for the parameter 
choice Ms = 10 13 GeV and vrl = 10 14 GeV. The mSugra parameters have been taken as in the 
SPS3 benchmark point and neutrino oscillation data have been fitted with the Y v fit, assuming 
degenerate right-handed neutrinos. 

note that, for very small seesaw scales all lepton flavour violating effects are negligible and 
no asymmetry is produced, since Al ~ Ar ~ 0. 

In addition, figure [7] shows again the relevance of Vr, which determines the parity breaking 
scale at which the LFV entries in the right-handed slepton sector essentially stop running. 
Lighter colours indicate larger Vr. As shown already in figure [6] for a particular point, the 
positron polarization approaches +1 as Vr is increased. 

Below the SU (2) r breaking scale parity is broken and left and right slepton soft masses 
evolve differently. The approximate solutions to the RGEs in equations ( 126|) and fl27|) show 
that, if neutrino data is fitted according to the Y v fit, the left-handed ones keep running 
from the SU(2)r breaking scale to the U(1)b-l scale. In this case one expects larger flavour 
violating effects in the left-handed slepton sector and a correlation with the ratio vbl/vr, 
which measures the difference between the breaking scales. This correlation, only present in 
the Y v fit, is shown in figure [HJ On the one hand, one finds that as vrl an d Vr become very 
different, Vbl/vr 1, the positron asymmetry approaches A = +1. On the other hand, 
when the two breaking scales are close, Vbl/ v r ~ lj this effect disappears and the positron 
polarization asymmetry approaches A = 0. Note that the Y v fit does not usually produce a 
negative value for A since the LFV terms in the right slepton sector never run more than the 
corresponding terms in the left-handed sector. The only possible exception to this general 
rule is, as discussed above, in the limit of very small m and Vrl/vr ~ 1. 

The determination of the ratio Vrl/vr from figure El is shown to be very inaccurate. 
This is due to the fact that other parameters, most importantly ttigut (which itself has an 
important dependence on the values of vbl an d Vr), have a strong impact on the results. 
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FIG. 7. Positron polarization asymmetry A(li + — > e + j) as a function of the seesaw scale, denned 
as the mass of the lightest right-handed neutrino, for the parameter choice vrl = 10 15 GeV and 
v R e [10 15 , 10 16 ] GeV. Lighter colours mean higher values of vr. The mSugra parameters have been 
taken as in the SPS3 benchmark point and neutrino oscillation data have been fitted with the Y v 
fit, assuming degenerate right-handed neutrinos. 
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FIG. 8. Positron polarization asymmetry A{^ + — > e + 7) as a function of the ratio vbl/vr- The 
seesaw scale Ms has been fixed to 10 13 GeV, whereas vrl an d vr take values in the ranges vrl G 
[10 14 , 10 15 ] GeV and vr G [10 15 ,10 16 ] GeV. Lighter colours indicate larger vrl- The mSugra 
parameters have been taken as in the SPS3 benchmark point and neutrino oscillation data have 
been fitted with the Y v fit, assuming degenerate right-handed neutrinos. 
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FIG. 9. Positron polarization asymmetry A(fi + — > e + 7) as a function of 

\og(v a/ 'trout) I ' ^og^BL/mGUT) ■ The parameters have been chosen as in figure[8l 
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FIG. 10. Positron polarization asymmetry A(/-t + — > e + 7) as a function of vr for three different 
mSugra benchmark points: SPSla' (blue line), SPS3 (green line) and SPS5 (red line). In this figure 
a fixed value vbl = 10 14 is taken. Neutrino oscillation data have been fitted with the / fit. 

Therefore, although it would be possible to constrain the high energy structure of the theory, 
a precise determination of the ratio vbl/vr will require additional input. Figure [91 on the 
other hand, shows that the polarization asymmetry A(fi + — > e +r y) is much better correlated 
with the quantity \og(vR/rriGUT) / log(t> blI^gvt)- This is as expected from equations ( 1261) 
and ( I27p and confirms the validity of this approximation. 

We close our discussion on the positron polarization asymmetry with some comments on 
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the / fit. Since this type of fit leads to Am 2 L ~ Arrive ~ in the Vrl — vr energy region, 
there is little dependence on these symmetry breaking scales. This is illustrated in figure [TUl 
where the asymmetry A is plotted as a function of Vr for three different mSugra benchmark 
points: SPSla' (blue line), SPS3 (green line) and SPS5 (red line). One clearly sees that the 
dependence on vr is quite weak compared to the Y v fit. In fact, the variations in this figure 
are mostly due to the changes in the low energy supersymmetric spectrum due to different 
vr values. In the case of the /-fit one then typically finds A G [0.0 — 0.3]. 

D. LFV at LHC/ILC 

Lepton flavour violation might show up at collider experiments as well. Although the 
following discussion is focused on the LHC discovery potential for LFV signatures, let us 
emphasize that a future linear collider will be able to determine the relevant observables 
with much higher precision. 

Figure [TT] shows Briji — > Xi e) and Br{fi — > x? A 4 ) as a function of the seesaw scale. The 
dashed lines correspond to T\ ~ tr and the solid ones to t^ — t^. As in the case of /i — > cy, 
see figure HI lower seesaw scales imply less flavour violating effects due to smaller Yukawa 
couplings. Moreover, figure [TT1 presents the same results for two different benchmark points, 
SPSla' and SPS3. As already shown in figure HJ /i — > is strongly dependent on the SUSY 
spectrum. For lighter supersymmetric particles, as in the benchmark point SPSla', fi — > e'-f 
is large, setting strong limits on the seesaw scale and thus on the possibility to observe LFV 
at colliders. In the case of heavier spectrums, as in SPS3, fi — > ej is still the most stringent 
constraint, but larger values of the seesaw scale and thus LFV violating branching ratios 
become allowed. Whether decays such as Br(fi — > Xi e ) an d Briji — > Xi A*) are observable 
at the LHC or not, thus depends very sensitively on the unknown m , M\ii and Mg. 

Furthermore, the right panel of figure [11] also shows that right staus can also have LFV 
decays with sizable rates. Of course, as emphasized already above, this is the main novelty 
of the LR model compared to pure seesaw models. This is direct consequence of parity 
restoration at high energies. 

Moreover, as in our analysis of the positron polarization asymmetry, one expects to find 
that if the difference between vr and vrl is increased, the difference between the LFV 
entries in the L and R sectors gets increased as well. This property of the Y v fit is shown in 
figure [T2| which shows branching ratios for the LFV decays of the staus as a function of Vrl 
for Vr G [10 15 , 5 • 10 15 ] GeV. As the figure shows, the theoretical expectation is confirmed 
numerically: the difference between Br(fi) and -Br(f^) strongly depends on the difference 
between vr and vrl- 

The question arises whether one can determine the ratio vrl/vr by measuring both 
Bt(tl) and Br{jR) at colliders. Figure [TBI attempts to answer this. Here the ratio Bt(tr — > 
Xi e)/Br(fi — > Xi e ) i s plotted as a function of vrl/vr. A measurement of both branching 
ratios would allow to constrain the ratio Vrl/vr an d increase our knowledge on the high 
energy regimes. For the sake of brevity we do not present here the analogous plots for 
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FIG. 11. Br(ji — > X\ e) and Br{ji — > Xi A 4 ) as a function of the seesaw scale, denned as the mass 
of the lightest right-handed neutrino, for the parameter choice vrl = 10 15 GeV and vr = 5 ■ 10 15 
GeV. The dashed lines correspond to T\ — tr and the solid ones to T2 — Tl- To the left, the mSugra 
parameters have been taken as in the SPSla' benchmark point, whereas to the right as in the SPS3 
benchmark point. In both figures neutrino oscillation data have been fitted according to the / fit, 
with non-degenerate right-handed neutrinos. The blue shaded regions are excluded by li — > e-y. 
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FIG. 12. Bt(tl — > X\ A*) an( i Br(fn — > X\ A 1 ) as a function of vbl for Ms = 10 13 GeV and 
vr € [10 15 ,5 • 10 15 ] GeV. Red dots correspond to T\ — tr, whereas the blue ones correspond to 
T 2 — T L- The mSugra parameters have been taken as in the SPS3 benchmark point and neutrino 
oscillation data have been fitted with the Y v fit, assuming degenerate right-handed neutrinos. 
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FIG. 13. Br(fji Xi ^)/Br{jL — > Xi A*) as a function of vbl/vr- The seesaw scale Mg has 
been fixed to 10 13 GeV, whereas vbl an d vr take values in the ranges vbl G [10 14 , 10 15 ] GeV and 
w/j £ [10 15 ,10 16 ] GeV. Lighter colours indicate larger vbl- The rest of the parameters have been 
chosen as in figure [T2l 
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FIG. 14. Br(f R X? e)/Br(f L Xi e) as a function of log(v R /m GUT )/ log(vBL/m GUT )- The 
parameters have been chosen as in figure [TH1 

other LFV slepton decays and/or other lepton final states, since they show very similar 
correlations with vrl/vr- For example, in principle, one could also use the ratio Br(fiR — > 
Xi r)/Br(jiL — > X? r ) to determine the ratio between the two high scales. 

However, as observed also in the polarization asymmetry for /i — > e-f there is an important 
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dependence on other parameters of the model, especially the exact value of rriGUT- This 
implies a theoretical uncertainty in the determination of vrl/vr. Again, as for A, a much 
better correlation with log(v r/itigut)/ log(i; bl/^gut) is found, see figure [HI 

In conclusion, to the determine vbl an d Vr individually more theoretical input is needed, 
such as the GUT scale thresholds, which are needed to fix the value of mcur- Recall, that we 
did not specify the exact values of these thresholds in our numerical calculation. This leads 
to a "floating" value of mcur when vr and vbl are varied. Also more experimental data is 
needed to make more definite predictions. Especially SUSY mass spectrum measurements, 
which may or may not be very precise at the LHC, depending on the SUSY point realized 
in nature, will be of great importance. Recall that, if in reach of a linear collider, slepton 
mass and branching ratio measurements can be highly precise. 

So far only slepton decays have been discussed. This served to illustrate the most inter- 
esting signatures of the model, namely, lepton flavour violation in the right slepton sector, 
however, LHC searches for lepton flavour violation usually concentrate on the decay chain 



75 



-77 



This well known signature has been widely studied due to the accurate information it can 



provide about the particle spectrum 78l-l82|. Note that in this decay one assumes usually that 
the X2 themselves stem from the decay chain qV, qx®- If the mass ordering m^o > mj > m^o 



is realized, the dilepton invariant mass 80|, [83|, defined as m 2 (l + l ) = (pi+ + Pi-) 2 , has an 



edge structure with a prominent kinematical endpoint at 

(m 2 — m 2 )(m 2 — m 2 ) 

where the masses of the charged leptons have been neglected. The position of this edge can 



be measured with impressively high precision at the LHC [781480 1 . implying also an accurate 
determination of the intermediate slepton masses. 

In fact, if two different sleptons / 12 have sufficiently high event rates for x\ ^iVJ ~~ * 
XilflJ an d their masses allow these chains to be on-shell, two different dilepton edge dis- 



tributions are expected |30l . l84j . This presents a powerful tool to measure slepton mass 
splittings, which in turn allows to discriminate between the standard mSugra expectation, 
with usually negligible mass splittings for the first two generations, and extended models 
with additional sources of flavour violation. 

The relation between the slepton mass splitting and the variation in the position of the 



2 2 — 4 

mi m- — m 



kinematical is edge is found to be |30l 

Amu Amj 



rriu mj (mj - m% )(mj - m 2 ) 

t A.1 t AO 



(36) 



Here Amu(i,j) = m^^ —m\ j i i is the difference between two edge positions, Amj = m~ v —m~ L 
the difference between slepton masses and fhu and m, average values of the corresponding 
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FIG. 15. ^Bi (left-hand side) and (right-hand side) as a function of the seesaw scale, defined 

win ni'i 

as the mass of the lightest right-handed neutrino, for the parameter choice vbl = 10 15 GeV and 
Vr 6 [10 15 , 10 16 ] GeV. Blue dots correspond to the mass distribution generated by intermediate left 
sleptons whereas red dots correspond to the mass distribution generated by the right ones. The 
mSugra parameters have been taken as in the SPS3 benchmark point and neutrino oscillation data 
have been fitted according to the Y v fit, with degenerate right-handed neutrinos. 



quantities. Note that higher order contributions of have been neglected in equation 



A number of studies about the dilepton mass distribution have been performed 



78-80 



concluding that the position of the edges can be measured at the LHC with an accuracy up 



to 10 . Moreover, as shown in reference |30|, this can be generally translated into a similar 
precision for the relative e — jl mass splitting, with some regions of parameter space where 
values as small as 10~ 4 might be measurable. Since this mass splitting is usually negligible 
in a pure mSugra scenario, it is regarded as an interesting signature of either lepton flavour 
violation or non-universality in the soft terms. Furthermore, in the context of this paper, it 
is important to emphasize that pure seesaw models can have this signature only in the left 
slepton sector (3l| . 

Figure [15] shows our results for the observables and as a function of the seesaw 
scale. Large values for Ms lead to sizable deviations from the mSugra expectation, with a 
distinctive multi-edge structure in the dilepton mass distribution. Moreover, this effect is 
found in both left- and right- mediated decays. Observing this affect would clearly point 
towards a non-minimal seesaw model, such as the LR model we discuss. 



As expected, these observables are correlated with other LFV signals 3jJ, |85j. Figure 



shows Br(fi ey) as a function of (^^f^J (mass distribution with intermediate L 
sleptons) and {^§^f^ (mass distribution with intermediate R sleptons). Again, the main 
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FIG. 16. Br(fj, —7- e^f) as a function of (^i^f") (left-hand side) and ^^Ft^) (right-hand side). 
The parameters are chosen as in figure [T5l 



novelty with respect to the usual seesaw implementations is the correlation in the right 



sector, not present in the minimal case 3l| • 

Furthermore, the process x\ ~~ * XiC^7 mi ght provide additional LFV signatures if the 
rate for decays with li ^ lj is sufficiently high. Reference (86| has investigated this possibility 
in great detail, performing a complete simulation of the CMS detector in the LHC for the 



decay x° — Xi^- The result is given in terms of the quantity 



K. 



Br(x° 2 -> X?e/i) 



eft 



Br(x° 2 -»■ X?ee) + Br(x° 2 XiW) ' 
which parametrizes the amount of flavour violation in X2 decays, 
the CMS test point LM1 (m = 60 GeV, M 1/2 = 250 GeV, A 
sign(fi) = - 



81 



(37) 

The study, focused on 
= GeV, tan/3 = 10, 



integrated luminosity of 10/6 



concludes that LFV can be discovered at the LHC at 5a level with an 
1 if K eix > K™ n = 0.04. 
Figure [T7] shows our computation of K efl as a function of the lightest right-handed neutrino 
mass, for the parameter choice vbl — 10 15 GeV and vr = 5 • 10 15 GeV. The results are 
shown splitting the contributions from intermediate left (blue) and right (red) sleptons. 
Although the selected mSugra parameters belong to thej3PS3 point, and not to LM1 as in 
reference 



86 



a similar sensitivity for K™ n is expected^. This is because the reduction in 
the cross-section due to the slightly heavier supersymmetric spectrum is possibly partially 
compensated by the corresponding reduction in the SM background and thus a limiting value 

uses 10 and with larger integrated 
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K™ n of a similar order is expected. Moreover, 
luminosities even smaller K™ n should become accessible at the LHC. 

The main result in figure [T7] is that for large Ms values the rates for LFV X2 decays 
are measurable for both left and right intermediate sleptons. In fact, for M5 > 10 12 GeV 



7 Moreover, the LM1 point, being very similar to SPSla', is strongly constrained by /i 
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FIG. 17. K e ^ as a function of the lightest right-handed neutrino mass, for the parameter choice 
vbl = 10 15 GeV and vr = 5-10 15 GeV. The blue curve corresponds to contributions from intermedi- 
ate L sleptons, whereas the red one corresponds to intermediate R sleptons. The mSugra parameters 
have been taken as in the SPS3 benchmark point, which satisfies mfy®) > m {h) > ^(x?)) an d thus 
the intermediate L and R sleptons can be produced on-shell. Neutrino oscillation data have been 
fitted according to the / fit, with non-degenerate right-handed neutrinos. The blue shaded region 
is excluded by fj, — > &y. 



the parameter K efJj is above its minimum value for the 5cr discovery of xi ~^ Xi e ^- See 



references 



this channel. 



86l . I87| for more details on the LHC discovery potential in the search for LFV in 



IV. CONCLUSIONS 

We have studied a supersymmetric left-right symmetric model. Our motivation for study- 
ing this setup was twofold. First, LR models are theoretically attractive, since they contain 
all the necessary ingredients to generate a seesaw mechanism, instead of adding it by hand 
as is so often done. And, second, in a setup where the SUSY LR is supplemented by flavour 
blind supersymmetry breaking boundary conditions, different from all pure seesaw setups, 
lepton flavour violation occurs in both, the left and the right slepton sectors. 

We have calculated possible low-energy signals of this SUSY LR model, using full 2-loop 
RGEs for all parameters. We have found that low-energy lepton flavour violating decays, 
such as fi — > e7 are (a) expected to be larger than for the corresponding mSugra points 
in parameter space of seesaw type- 1 models and (b) the polarization asymmetry A of the 
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outgoing positron is found to differ significantly from the pure seesaw prediction of A = +1 
in large regions of parameter space. We have also discussed possible collider signatures of the 
SUSY LR model for LHC and a possible ILC Mass splittings between smuons and selectrons 
and LFV violating slepton decays should occur in both the left and the right slepton sector, 
again different from the pure seesaw expectations. 

We think therefore that the SUSY LR model is a good example of a "beyond" minimal, 
pure seesaw and offers many interesting novelties. For example, the impact of the interme- 
diate scales on dark matter relic density and on certain mass combinations and the influence 
of the right-handed neutrino spectrum on low energy observables, are topics that certainly 
deserve further studies. 
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Appendix A: RGEs 

We present in the following appendices our results for the RGEs of the model above the 
U(1)b-l breaking scale. We will only show the /3-functions for the gauge couplings and the 
anomalous dimensions of all chiral superfields. We briefly discuss in this section how these 
results were calculated. Furthermore, we show how they can be used to calculate the other 
/3-functions of the models and give as example the 1-loop results for the soft SUSY breaking 
masses of the sleptons. The complete results are given online on this site 

http : / / theorie . physik . uni-wuerzburg . de/~f nstaub/supplementary . html 

In addition, the corresponding model files for SARAH are also given on this web page. 
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1. Calculation of supersymmetric RGEs 



For a general N = 1 supersymmetric gauge theory with superpotential 



W(<f>) 



2 r~ t.tj g 

the soft SUSY-breaking scalar terms are given by 



soft 



1 



b x3 <t>4j + -h l3k ^^ k + c.c. ) + (myjfacj)* . 



(Al) 



(A2) 



The anomalous dimensions are given by 



51 



7 



(i)j 



7, 



(2)j 



1 



Y imn Y n ^Y pqr Y mr i + ( 7 2 F w y^^[2C , 2 (p) - C 2 (i)} 



(A3) 



(A4) 



+ 25^ 4 [C 2 (z)5( J R) + 2C 2 (z) 2 - 3C 2 (G)C 2 (z)] , 

and the /3-functions for the gauge couplings are given by 

/?« =g 3 [S(R) - 3C 2 (G)} , (A5) 
=9 5 {-6[C 2 (G)} 2 + 2C 2 (G)S(R) + 4S(R)C 2 (R)} - g 3 Y l3 % k C 2 (k)/d(G) . (A6) 

Here, C 2 (i) is the quadratic Casimir for a specific superfield and C 2 (R),C 2 (G) are the 
quadratic Casimirs for the matter and adjoint representations, respectively. d(G) is the 
dimension of the adjoint representation. 

The /3-functions for the superpotential parameters can be obtained by using superfield tech- 



nique. The obtained expressions are 88|, [89 



Bf = Y^jV , 



'Y 

Qij 

7' 



7 P 



(AT) 
(A8) 



The (..) in the superscripts denote symmetrization. Most of the B-f unctions of the models 
can be derived from these results using the procedure given in [90] based on the spurion 
formalism 
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In the following, we briefly summarize the basic ideas of this calculation for 



completeness. 



The exact results for the soft /3-functions are given by 90] : 



M 



20 



09 

9 



{B m ,y 3 = ^ 3 . 



Bl j = b l ^\ - , 



(A9) 



(A10) 
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where we defined 



O = Mg 2 — - h lmn —— , (All) 
(71ft = 1 i j , (A12) 



d 

A = 200* + 2MM*g 2 



dg'< 



Y lmn -^— + c.c. 



+ X§-. (A13) 



the known contributions of X from 
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Here, M is the gaugino mass and Y ljk = (m 2 ) l iY jkl + (m 2 y{Y lkl + (m 2 ) k iY ijl . Eqs. flA9[) - 
(1A10I) hold in a class of renormalization schemes that includes the DRED'-one (92|. We take 



X DRED'(1) = _ 2g 3 S ? (A14) 

x dred'( 2 ) = (2 r )- x ghr[WC 2 {R)\ - Ag 5 C 2 (G)S - 2g 5 C 2 (G)QMM* , (A15) 

where 

S = r~hi[m 2 C 2 (R)} - MM*C 2 {G) , (A16) 
W\ = ^Y ipg Y^ n (m 2 y n + ]Y^Y vqn {m 2 Y i + 2Y ipq Y^{m 2 )\ + h ipq V™ - 8g 2 MM*C 2 (Ry i . 

(A17) 

With Q = T(R) — 3C 2 (G), and T(R) = tr [C 2 (R)], r being the number of group generators. 

2. Prom GUT scale to SU(2)r breaking scale 
In the following sections we will use the definitions 

Yq{ = , Yii=yi ik (Ais) 

and in the same way Tq and T l £ . We will also assume summation of repeated indices. 
a. Anomalous Dimensions 



^ = 2Y Q k Y S k - y 2 (l8</ 2 2 + 32^ 2 + 9bl) 1 (A19) 

7 } = + m ( " 128 ^3 + 2052 ^ 4 + 289 ^ + 3% 2 ( 32 ^ + &l) + ^g 2 g 2 BL ) 1 
+ y* Qm (Zg 2 2 5 mn + ^Tr (aa*)6 mn - 2Tr (^lf m ) - 6Tr (y^ Y^))yI 



32Yg m Yq n Yq k Yg m (A20) 
2Y k Y Q* ~ J2 i 1892 + 32 ^ + ^ 1 (A21) 
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7 (2) 



7 



^- ( - 12% 3 4 + 2052^ + 289<?| L + 36g 2 2 (32g 2 3 + «?| L ) + 64^I L ) 1 

+ (e A» + - 2Tr (V;^) " 6^^^))^ 

- 32 >Q m ^^Q r m ^ (A22) 
(1) " 2 (3/t/ + YtYl) - \ (2gj + g% L ) 1 (A23) 



1l 16 



12^| L + 76<? 2 4 + 99^ L )l + 3/ t /( - 3|a| 2 - 4Tr(ff) + 6 ^ + 8 ^ 
+ ^ (6 An + ^Tr (aa*)^ - 2Tt(y*Y^) ~ 6Tr (y^ 



n/V^n)^; - 4y; m yi n y in yl - 2/t(i?//t + w Lk vl)f 

S?YlJYL (A24) 



(i) 



2(3/^ + ^^^) - i (2^ 2 2 + ^ L )l (A25) 
(\2glg 2 BL + 76^ 2 4 + 99^) 1 + 3// t ( - 3|a| 2 - 4 r Lr(// t ) + 6g 2 BL + 8<? 2 2 ) 



+ Yl (Q9 2 2 S mn + -Tr(aa*)6 mn - 2Tr(Y*Y^ m ) - QTt(y*Y^ 

- llff5 mn )Y Ln - 4Yt m Y Ln Y[ m Y£ n - 2/(l7/V + SY^Y^f 
-QYlY L jf (A26) 

(7?^ = -3S 2 2 1 - §(aa' + a*a) + S im 5 jn (sTr (y$J&) + Tr (y^Y*)) (A27) 

(74 = 33g 4 l — 9^2(aaa*a* + «*a*aa) + 3(aa*aa* + a*aa*a)) 

- 24(aa* + a* a) (2$ + 2Tr ( 



T 



- 1 (*jm< + «> m ) (3Tv(y4 Yl) + Tr (y; m y L , 

- ^<W^( - 3g 2 BL Tv{Y* Lm Yl) - (32gj + gl^^Y^ 
+ 2(51v(//tY L ^ ^ 

+ 4(2Tr(y L t m y in yf n y; n ) + Tr (y^Y^) +^{YlY Ln YlYl n 
+ 12 (2Tr(y^ Y Qn YlY* n ) + Tr (y^Y^ Y?Y* m ) 

+ ^(yL Y Q^ Y L Y L))) (A28) 
7 « = 2Tr(//t) - 3g 2 BL - Ag 2 2 + \\a\ 2 (A29) 
7 i 2) = 48^ + 24^ + 81^ 

+ ^|a| 2 (4s 2 + Tr(a«*) - I|a| 2 ) - (2^ + 3g BL ) Tr (ff 
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- 24Tr(//t//+) - 6Tr(//ty Lfc y L t fc ) _ 2 Tr (/F^F^) (A30) 
^ = -^Il-^U\h 2 (A31) 
lf = |(4(l6^ + 27^ + 8^ £ ) 

+ |a| 2 (2Tr(cta*) -STr^) -7|a| 2 + 8# 2 )) (A32) 

7^ = 2Tr(//t) _ 3^ _ 4 ^2 + | |a|2 {m 

7 g = 48<? 2 4 + 24<&| L + Slg% L 

+ ^H 2 (4^ 2 + Tr(aa*) - ^|a| 2 ) - (2g\ + S^Tr (//t) 

- 24Tr(//t // tj _ 8Tr(/lf fc F i : fc /tj (A34) 
7i 1 c ) = -3^ L -4^ 2 + ^|a| 2 (A35) 
72 = ^(4(16^ + 27^ + 8^1,) 

+ |a| 2 (2Tr(aa*) -3Tr(// f ) -7|a| 2 + 8^)) (A36) 

7 j» = 2\a\ a - 4g* (A37) 



(2) 



(2) 



3Tr(a(aa* - + 48^ + |a| 2 (l2^ L - 3Tr (// f ) -6|a| 2 + 8#f) (A38) 

2\a\ 2 -Ag 2 2 (A39) 
3Tr(a(aa* - + 48^ + |a| 2 (l2^ L - 3Tr (// f ) -6|a| 2 + 8#f) (A40) 



Note that the previous formulas are totally general and can be applied with any num- 
ber of bidoublets. Nevertheless, if two bidoublets are considered aa* = a* a and further 
simplifications are possible. 



b. Beta functions for soft breaking masses of sleptons 

Using the procedure explained in sec. IA 11 we can calculate the soft breaking masses for 
the sleptons. The results are 
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where 



IQn^ml = 6ffW L + UfmlP + Qm 2 L ff + \2m\f f 



dt 



-2Y Lk Ylm 2 L + 2m\Y Lk Yl + 4Y Lk m 2 Lc Yl 



+±(ml) mn Yt ) Yt ) ' + UT f T} + AT Lk Tl 
-(S^jMil 2 + 6gj\M 2 \ 2 + \g 2 BL S x )l (A41) 
l 6vr 2^ m 2 c = 6/ t /m 2 c + i2/t m | c/ + 6m| c /t/ + 12ml e pf 



dt 



+2YlY Lh m 2 LC + 2m 2 LC YlY Lk + ^lra 2 L Y Lk 
+^ml) mn Yt ] V W + 127*7) + 4Tt T ifc 

"(S^l^il 2 + 6^ 2 |M 2 | 2 - 5^)1 (A42) 



Si = 3(m A — — m A c + m 



+ [K)mn - (™Qc) mn ~ (™Dmn + (^L)mn] (A43) 



c. 5eia functions for gauge couplings 



til = ^9 3 BL (A44) 

til = ~ 192|a| 2 - 93Tr(//t) + 2(ll5^| L + 162^ 

" 2Tr(l^ y^) - 6Tr(V L * r£) + 8 ^ 2 )) (A45) 

4 1] = 8sf (A46) 
= ^2 3 (660^ 2 2 + 144^ + 162^| L - 192|a| 2 + 108Tr(aa*) - 73Tr(//t) 

- 24Tr(F£ fc y£) - 72Tr(^ fc l^)) (A47) 

= -%3 (A48) 

P$ = gl (Ug 2 3 + 18g 2 2 - 8Tr (y^) + g 2 BL ) (A49) 



3G 



3. Prom SU(2)r breaking scale to U(1)b-l breaking scale 



a. Anomalous Dimensions 

- 2(r;yJy;rJ + y u *yJY u *y w T ) 

+ W T ( - 3Tr(F d Fi) - |6 C | 2 - \\b\ 2 - Tr(F e F e t) + , 2 ) 
+ Y:YZ( - 3Tr(F u F u t) - \b c \ 2 - ^| 2 - Tr(r^t) + ^) 
T « = 2y;y d - 1 (32^ 2 + 6^ 2 , + 1 

^ = - lii ( - 109 ^ + 128 ^ 4 - 12 9bl9 2 r ~ 64^| (6^ 2 + g 2 BL ) - 684^) 

- YjY d (2\b c \ 2 + 2Tr(Y e Y^ + 3|6| 2 - Qg 2 L + 6Tr(F d y;)) 
7 W=2F u tF tt -l(32, 2 + 6, 2 +,| L )l 

^ = ~ 144 ( " 109 ^ L + 128 ^ " 12 ^ L ^ " 64 ^ ( 6 ^ + " 684 ^«) 

-2(yJ-y d y]y u + y u ty u y u ty u ) 

- y,Jy u (2|6 c | 2 + 2Tr(Y v Yj) + 3\b\ 2 - 6g 2 L + QTt(y u Y^)) 

if = -\ (2yi + </L) i + n*n T + n% T 

7 f = + A (124 + i3^ L + 4s 2 L g\ ) i - 2 (y;y/y;y e T + y;yJX% T ) 

- y;f c y u t 



+ y e *Y e T ( 


^-3Tr(y d y;)- 


N 2 - 


>l 2 


-Tr(y e y e t) 




+ Y u *Yj( 


^-3Tr(y u y u t)- 




>' 2 






= 2YjY e 


- J (2<7 2 + 3<^ 


)1 









7$ = (ll7</| L + 12^ L ^ 2 + 76^) 1 " 2 (^^^ + ^ Y V Y*Y^ 

- y e t y(2|6 c | 2 + 2Tr(yy e t) + 3|&| 2 - 6 9 2 L + 6Tr(y d y;)) 

7? = 2YX - \ (2^ + 3g 2 BL ) 1 + F c * 
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- 2yj(y e y e t + yyj)y - 2(/ c 1 + + f l c *)Y?Y:ul + / C 1T ) 

- y;y„(2( - 3 5 | + 3Tr(y u y u t) + |& c | 2 + ^(y.yj)) + 3|6| 2 ) 

+ ^ (2^ + 3g 2 BL - - Tr(fj?) - Tr^t/f)) (A61) 
T g = 3Tr(y d y;) - 1^ + 3 -\bf - \g\ + |6 C | 2 + Tr(y e I?) (A62) 

+ |6| 2 ( - 18Tr(yy u t) - 15|6| 2 - 19|6 C | 2 + 24g 2 L - 6Tr(yyj)) 

- 4|6 C | 2 (3|6 C | 2 + 3Tr(yy u t) + \a\\ 2 + 1i(y u Y})) ~ 36Tr (Y d YjY d Y^ 

- 12Tr (Y d YjY u Y^ - 12Tr (V^y^) - 4Tr (yy^yj) ) (A63) 
7 « = 3Tr(yy u t) - 1^ + ^|6| 2 - ^1 + |6 C | 2 + Tr(yyj) (A64) 

^S = i( 27 ^ +6 ^ +19 ^ 

+ |6| 2 ( - 18Tr(yy;) - 15|6| 2 - 19|6 C | 2 + 24^ - 6Tr(yy e t )) 

- 4|6 C | 2 (3|6 C | 2 + 3Tr(yyj) + \a\\ 2 + Tr(y.l?)) + 2(32^ + & l )Tc(yX) 

+ 6^| L Tr(yy 1 t) -^(/i/ityjy;) -4Tr(/ c 1 y;y/ c 1 t) _ i 2 iy (y d y]y u y u t) 

- 4Tr(y e y e tyy;) - 36Tr (y u yjy u y u t) - 4Tv(y/ c 1 t/f yj) 

- i2Tr (yyjyyj) - 4iv(/ c 1 tyjy;/f)) (A65) 

T^o = -2^ - 3<£ x + K I 2 + Tr + Tr (/^/f ) (A66) 

7go = ~ (72^x + 24^4 + 444 - A\al\ 2 (|a^| 2 + |6 C | 2 ) 

- (24 + 3,L) (Tr^/^t) + Tr^t/f )) - ATr(f c fJf c fj) 

-x*(fif?fl; T ff) -™{flf?fl T fl*)-^{flfMY:) 

- ATr^YXfj) - 4Tr(y/ c 1 t/f yj) - 4tt(#t/f /*t/f ) 



- 8Tr(/^t/f ^/f) - ^(^tyjy;/^ )) (A67) 
7 «=-24 -3sI L + |a*| 2 (A68) 
7? i = 12<7L4 + 224 + 36g 4 BL 

- |^| 2 (2|^| 2 + 2|6 C | 2 -h Tr(/ c Vit) + ^(/^t^^ (A6g) 

7? = -4^1 + |&| 2 (A70) 
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7 < 2) = 28g 4 L - \b\ 2 (2\b c \ 2 + 3\b\ 2 + 3Tr(Vj) + STrjV^) 

-gl + gl + Tr (y e y e t) + Tr (Y u Yj) ) (A71) 
l£ = 2\b c \ 2 + K\ 2 (A72) 
it = -2|^| 4 - |aj| 2 ( - 4<& - 6,| L + Tr^/^t) + I^t/f )) 

- 2|6 C | 2 (2|6 C | 2 + 3|6| 2 - + 3Tr(Y d yj) + 3Tr(V u Y^ 

- ^ + Tr (y e y e t) + Tr (y v YJ) ) ( A73) 

In these expressions we have defined 

F c = flfl ] + fife* + fl T fl ] + /f/c 1 * (A74) 

b. Beta functions for soft breaking masses of sleptons 

Again, the results for the slepton soft SUSY breaking masses at 1-loop are shown. The 
beta functions read 

lQK 2 jmi = 2F e m 2 c y e t + 2m 2 H Y e Y^ + 2m 2 H Y v Y^ + m 2 L Y e Y^ 
+Y e Y^m 2 L + m 2 L Y v Yl + Y u Y^m 2 L + 2Y v m 2 v0 Y^ 
+2T e Tl + 2T v Tl - (S^lJMxl 2 + 6g 2 L \M L \ 2 + ~g 2 BL S 2 )l (A75) 

16^. = 2YXm% + 2m 2 ^ e Y e + 4<>? y. + ^m 2 L Y e 

+ATlT e - (3g> BL \Mi\ 2 + 2g 2 R \M R \ 2 - ^g 2 BL S 2 - \g 2 R S 3 )l (A76) 

where 

5 2 = 2{m\ c o - m\ c0 ) + Tr [m 2 dc - m 2 c + m 2 c - m 2 c + 2m 2 - 2mg] (A77) 

5 3 = 2{m\ c0 - m\ c0 - m 2 Hd + m 2 H J + Tr [3m^ c + m 2 c - 3m 2 c - m 2 c ] (A78) 

c. Beta functions for gauge couplings 



= $9 3 bl (A79) 
P£L = + W9 2 bl + Wl + 30^1 - 12|a^| 2 - QTr^ 1 /^) - 2Tr(yyj) 

- 6Tr(r e F e t) - 2Tr(y u yJ-) - 6Tr(yy,t) - gTr^Vf )) (A80) 

P$> = 3^1 (A81) 
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9L 



PS 



PS 



i 9 2( - 28|6| 2 + 3(24 fc 2 - 2\bf - 21t(y,^ - 27r(\'X) + 3» 

6Tr(y i F])-6iv(y'»y;t)+ J 2 i 



2 

BL 



90* 

g 3 R (24g 2 3 + 15g 2 BL + 3g\ + 15g 2 R - 4|^| 2 - 4|6| 2 - 2|6, 



3Tr 



6iv(y d y d t ) - 2Tr(r e F e t) - 6Tr(Y u r u t) _ 2Tr(ni;t) - 3Tr^ff 
-3^ 3 3 

<7 3 3 (l4^ 3 2 + 3g 2 R - 4Tr (f,^) - 4Tr (y u Y^) + 9g 2 L + g 2 BL ) 



(A82) 
(A83) 

(A84) 
(A85) 
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